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Abstract 

We study the QCD evolution for the twist-three quark-gluon correlation functions associated 
with the transverse momentum odd quark distributions. Different from that for the leading twist 
quark distributions, these evolution equations involve more general twist-three functions beyond 
the correlation functions themselves. They provide important information on nucleon structure, 
and can be studied in the semi-inclusive hadron production in deep inelastic scattering and Drell- 
Yan lepton pair production in pp scattering process. 
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Semi-inclusive hadronic processes have attracted much theoretical interest in recent years, 
where the so-called transverse momentum dependent (TMD) parton distributions and frag- 
mentation functions can be studied naaaaia , 13, |8j. These functions generalize the 

original Feynman parton picture, where the partons only carry longitudinal momentum 
fraction of the parenting (final state) hadron. They certainly provide further important 
information on nucleon structure, and are crucial to understand the novel spin phenomena, 



such as the single transverse spin asymmetry (SSA) [l|, |2|, |3|, |9|, 1 1 CM ] - 

Important aspects of the TMD parton distributions have been explored in the last few 
years, such as the gauge property and the crucial role of the initial/final state interaction for 
the nonzero Sivers quark distribution leading to the SSA in Semi-inclusive hadron production 
in deep inelastic scattering (SIDIS) and Drell-Yan lepton pair production processes 

. Further study has shown that this mechanism is uniquely related to the twist-three 

151 ] . In particular, 
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12 



13, 



14. 



quark-gluon correlation approach for the SSA phenomena 
these two approaches are unified to describe the same physics in the overlap region where 
both apply |l6 |. 

At the leading order, there are eight independent TMD quark distributions, depending 
on the polarizations of the nucleon and the quark [?], Three of them are called /el- 
even distributions. After integrating over transverse momentum, they produce the leading- 
twist quark distributions, including the spin average, longitudinal polarized, and transversity 
quark distributions [TZj. The rest five distributions are called fcj_-odd distributions. Upon 
integral over the transverse momentum, they will vanish in the quark correlation matrix. In 
this paper, we are interested in four of these fcj_-odd TMD quark distributions. They can 
be defined from the following matrix, 



M at3 (x,k A 



2tt 



(2*) 



PS 



a£(tw(r,o,&j.) PS ) » (!) 



where x is the longitudinal momentum fraction of the proton carried by the quark and k± is 
the transverse momentum. In the above definition we have chosen P = (P + ,CT,0j_) which 
is along the momentum direction of the proton, S is the polarization vector, and & v [£) is 
defined as & v (0 = C v (— oo; , with C v the gauge link. This gauge link contains the 

light-cone gauge link contribution and the transverse gauge link contribution at the spatial 
infinity [sj]. We have chosen it goes to — oo, indicating that we adopt the definition for the 
TMD quark distributions for the Drell-Yan process [l[ [2, [J. The four fc^-odd TMD quark 
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distributions can be obtained by the following expansion of the above quark correlation 



matrix 



SB, 



M 



1 



9it(x, k ± ) lb P(k ± ■ S ± ) + tf T (x, k A _y™P 1 »P u k a S l 



2M 

+h 1L (x, k^Xia^P^kl + hi(x, k ± )a^k^P u + •••], (2) 



where M is the nucleon mass and the interpretations of the four fc^-odd TMD quark distribu- 
tions are: g%T and represent a longitudinal polarized and unpolarized quark distributions 
in a transversely polarized nucleon, respectively; tin and represent transversely polarized 
quark distributions in a longitudinal polarized and unpolarized nucleon target, respectively 
They are /cj_-odd distributions, i.e., after integrating over the transverse momentum, the 
above expansion will vanish. However, if we weight the integral with transverse momentum, 
the above matrix will lead to a set of quark-gluon correlation functions at the twist-three 
level. These correlation functions can be calculated as transverse momentum-moment of the 
above four k±-odd TMD quark distributions. The last k±-odd TMD quark distribution hj~ T 
represents a correlated transversely polarized quark distribution in a transversely polarized 
nucleon target, and is related to the twist-four quark-gluon correlation function. We will 
not discuss this function in this paper. 

As mentioned above, the transverse-momentum-moment of the above four TMD quark 
distributions define the following transverse momentum dependent correlation functions in 
nucleon, 

k 2 f k 2 

d 2 k±^^f^ T {x,k ± ) = T F (x) , / d 2 k ±1 ^gi T (x,k ± ) = g{x) , (3) 

We emphasize again that the above TMD quark distributions follow their definitions in 
Drell-Yan process. If we choose those for the semi-inclusive DIS process, the above two 
equations associated with Tp{x) and Tp\x) shall change signs. These correlation functions 
are related to more general quark-gluon correlation functions. For example, Tp(x) and 
Tp\x) are diagonal parts of the general quark-gluon correlation functions Tj?(x-\ , X?,) and 



Tp(x±, £2) which are responsible for the single spin asymmetry in hadronic process 12, 18]: 
T F (x) = T F (x,x) and Tp(x) = Tp r \x,x). 1 They can be defined through the following 



1 For the convenience of our presentation, we have changed the normalizations for Tp{x\,X2) and 
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correlation matrix, 



M ^ x) = i d £ d -t elxP+y ~ 



(PS\MO)gF + ^y^ a (y-)\PS) 



(5) 



where \i is a transverse index, F +IJL the gluon field tensor, and the gauge link has been 



suppressed. Its decomposition contains the contribution from Tp(x) and T J 



X 



M 



T F {x)e^S^ + T { ;\x)i^J 
Similarly, we can calculate the other two correlation functions by {4] 



(6) 



M 

T 



(7) 



where Mp a g is defined as, 



2vr 



i£-xP+ 



dCgF + »(C)}MO\PS} (8) 



Applying the time-reversal invariance, we find the above definition of g is the same as that 
in [14j, except a normalization factor 2. 

The above four correlation functions are subsets of more general twist-three quark-gluon 
correlation functions [l^, 3]: Gn(x,y), Go(x,y), H D (x,y) and E(x,y). These twist-three 
functions and their contributions to the inclusive DIS and Drell-Yan lepton pair productions 
have been under intense investigations in the last two decades (see for example 17|). The 
above four correlation functions Eqs. (3,4), however, will enter in the transverse momentum 



weighted cross sections in the semi-inc 



usive hadron production in DIS and Drell-Yan lepton 



pair production in pp collisions [3, Is, H]- They will provide additional information on 
the quark-gluon correlations in nucleon, and will be complementary to those studied in 
the inclusive DIS and Drell-Yan processes. Recent experimental developments will help to 
pin down these contributions, and build strong physics associated with these correlation 
functions 21]. 



One of the important questions remained to be answered is the scale evolution for these 
correlation functions. The evolution equation controls the energy dependence of the associ- 
ated observables 22|. For example, with the evolution equations, we will be able to compare 
the single spin asymmetries coming from the same quark-gluon correlation function Tp{x) 



Tp (x\,X2) by a factor of 1/2-ttM as compared to those in |l6l. Il8|. 
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in hadronic processes at different energy experiments. General structure of the evolution 
equations for the twist-three quark-gluon correlation functions has been known in the liter- 
ature 23[]. However, the above correlation functions Eqs. (3,4) are special projections of the 



general twist-three quark-gluon correlations, and their evolutions are not directly available 
from the already known results 23]. Earlier attempts 24J] have been made to derive the 
evolution equations for the correlation functions of Eqs. (3,4), but were not complete. On 
;he other hand, from the large transverse momentum quark Sivers function calculated in 
1(|, we would already obtain the evolution equation for Tp(x), since the collinear diver- 
gence in that calculation will lead to the splitting function of Tp(x). This splitting function 
was confirmed by a complete calculation of next-to-leading order QCD correction to the 



transverse-momentum weighted spin asymmetry in Drell-Yan lepton pair production [25]. 
More comprehensive evolution equations for Tp(x), together with those for the three-gluon 
correlation functions which are relevant to the single spin asymmetry observables have re- 
cently been derived in 26| . In this paper, we will extend these studies to calculate the scale 
evolutions for the above four quark-gluon correlation functions. Their contributions to the 
azimuthal angle distributions in Drell-Yan lepton pair production in pp collisions, and the 
relevant QCD factorization analysis will be presented in a forthcoming publication 

In our calculations, we will choose the light-cone gauge: A + = 0. There are several 
advantages for this choice. First, the quark Sivers function was previously calculated in 



the covariant gauge la]. Our calculation in the light-cone gauge will provide an important 



cross check for the results. Second, the light-cone gauge is more convenient to calculate 
the evolution equations for g and h. In particular, the evolution equations for Tp(x) and 
g(x) can be calculated simultaneously. The only difference is that for Tp(x) we have to 
take a pole contribution for some diagrams, whereas for g(x) we will not take the pole (see 
the discussions below). Third, we can further choose a particular boundary condition in the 
light-cone gauge [3]] , which will greatly simplify the derivation. We have also checked that the 
final results do not depend on the boundary condition. According to the quark distribution 
definition we have chosen above, it is convenient to choose the retarded boundary condition, 
i.e., A±(— oo~) = 0. With this choice, the gauge link associated with the TMD quark 
distributions in Eq. (1) becomes unit, and their contributions can be neglected [3J. From 
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FIG. 1: Real gluon radiation contribution to the evolution equations for the twist-three quark-gluon 
correlation functions Tp(x), g(x), Tp (x), h(x). 



this, we can re- write the quark-gluon correlation functions Tp(x) and g as 



T F (x) 



9{x) 



dy 



8ir 2 M 



(PS$(0)ieTS ±u id ± ^ a (y-)\PS) 



dy 

4vr M l 



jixP+y 



{ps\Mo) l5 is ± ^My~)\PS) , 



(9) 
(10) 



in the light-cone gauge with retarded boundary condition. Similar expressions hold for 
other two correlation functions, and h. In the following calculations, we will focus on 
the derivation for the evolution functions for Tp and g, especially for Tp, and those for T F °^ 
and h can be obtained accordingly. 

To calculate the splitting function for the above two functions, we have to take into 
account the contributions from the operators {^pd±ip) and ^A±ijj^ [l9], because they are 
at the same order. Especially, because of the contribution from A±, the evolution of the 
above correlation functions will involve more general twist-three functions: Go and Gp> 
or T F (xi,x 2 ) and T F cr \xi,x 2 ). This is an important feature for the scale evolution of the 



23|]. 



higher-twist distributions, such as that of the gp structure function 

We plot the Feynman diagram contributions from the real gluon radiations in Fig. 1, 
where (a) is the contribution from the partial derivative on the quark field, and (b — d) are 
those from A± contributions. The virtual corrections only contribute to partial derivative 
part, and they are easy to carry out. 

We will perform the collinear expansion for the hard scattering part to calculate the 



contribution from Fig. 1(a). The linear k± expansion term combining with the quark field 
will lead to the quark-gluon correlation function T F (x) and g(x) in Eqs. (9,10). In the 
collinear expansion in terms of k±, we can fix the transverse momentum of the probing quark 
(Iq) or the radiated gluon (l g ), because of momentum conservation and we are integrating 
over them to obtain T F (x) and g(x). We have also checked that they will generate the same 
result. In the following calculations, we choose l g being fixed in the collinear expansion. 
This will avoid the collinear expansion of the on-shell condition for the radiated gluon, and 
simplify the derivations. 

For the A± contribution, we notice that F +M = d + in the light cone gauge. Therefore, 
one can relate the corresponding soft matrix to the correlation function Tp{x,X\) in the 
following way, 



&!) dUl c ixiP+y- e i(x- Xl )P+y- 



x — X\ + it J Att 
dy-dy 1 



{PS\MO-)ie^S ± „gF + M)My-)\PS) 



^ P+ e -^ + y- e *(--OP + yr(p^|^( -)^ ± ^ ±M (yr)^(y-)|P5) . (11) 

In the above formula, the soft gluon pole appears in the first line comes from the partial 
integration. The prescription of this pole has been determined because we have chosen 
the retarded boundary condition. For the same reason, we have to regulate the light cone 
propagator in a consistent manner, and the gluon propagator in Fig. 1(c) in the light cone 
gauge with retarded boundary condition is given by 3j, 

d*w - ^ (** - , (12) 

I 2 + it \ I ■ n + ie J 

where I is the gluon propagator momentum entering the quark-gluon vertex in Fig. 1(c). 

Adding the contributions from the partial derivative and A±, we reach our final formula 
for the Tp(x) splitting calculation, 

T F l \x B ) = j dxd% ± ^ {[H(k,l g )fl x \ k±=0 T F (x,x) 

+ / dxd Xl d% ± {[H tl {xP,x 1 P,l g )^>] x ZjJ - _ 1 T F (x, Xl ) , (13) 

where the transverse spin vector has been integrated out, and the transverse index \x is not 
meant to be summed up. l q ± is the probing quark transverse momentum. In the above 
equation, H(k,l g ) represents the hard partonic part in Fig. 1(a) with transverse momen- 
tum dependence on k±, and H^(xP,XiP,l g ) the hard part for Figs, l(b-d) with transverse 



polarized gluon insertion where all momenta are collinear. We have to include both 
contributions to obtain a complete result. 

We have similar expression for g(x) splitting. The only difference is to replace^ with 75 fa 
and similar replacement in the hard parts in the above two terms. More over, because of 
simple Dirac algebra, the first term is the same for both Tp(x) and g(x) which comes from 
Fig. 1(a). Let us first discuss this contribution. In the calculations, we have to perform the 
collinear expansion in terms of k^_. Because of the momentum conservation, l£ = k^_ — l£, 
we can separate the contribution from the explicit dependence on kj_, 

Tjp\mg.i(*) = J dxT F (x,x)d 2 l g± j H(k,l g )ft \k ± =o ~ l g±^F H(k,l g )$ | fc± =o| -(14) 

The first term is easy to derive, and its contribution will be 

a s fdxdl 2 g± f i + z ^ 

where z = xb/x, and the well-known splitting kernel appears. This splitting kernel contains 
the end-point divergence, which should be canceled out by the virtual diagram contributions. 
After taking into account the virtual contribution, the end-point will be regulated by the 
plus function, 



Oi, 



<h ' dl kc F ( ^±^- + - z) ) T F (x, x) , (16) 



2vri x P g± \(l-z) + 2 



where the plus function follows the definition of [22|. To calculate the second term of 
Eq. (14), we will do the collinear expansion of the hard scattering part H(k,l g ) with l g j_ 
fixed. The transverse momentum k± flow can go through the quark line in Fig. 1(a), for 
which we label as "cross" in the diagram. We can further simplify the derivation by using 
the following identity, 

6 1 %< , (17) 



dk a ^ ft ' It 

which essentially represents the application of the Ward identity. Applying the above iden- 
tity, we can relate the k± expansion in the quark propagator and quark line to that with a 
transverse polarized gluon insertion with zero momentum attachment. These contributions 
shall be combined with those from Figs, l(c-d). We will discuss them below. 

As we mentioned above, the contributions from Fig. 1(a) are the same for the evolutions 
of Tp{x) and g(x). Therefore, the above results apply for that of g(x) too. However, the 
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contributions from Figs, l(b-d) are different for T F (x) and g(x). This is because for this 
part, we have to take pole contribution to obtain the splitting for T F (x), whereas for g(x) 
we do not take pole contribution. We first discuss their contributions to the evolution of 



T F (x) correlation function. Depending on the value of x g 



X\ when we take the pole, 



these poles are called soft [x g = 0) and hard (x g 7^ 0) poles, respectively. The hard pole 
contribution only comes from the light-cone propagator in Fig. 1(c), and its contribution is 
easy to calculate. For this, we obtain 



T; 



(l)|hp 



F iFig.l(c) 



a s 
2ti 



dx dl 2 g± C A fl + z 



x P g± 2 



T F (xz, x) 



(18) 



We emphasize that for the hard pole contribution the explicit factor l/(x — x±) has been 
included in the above result, which is finite because x ^ X\. 

On the other hand, the soft pole contribution comes from the explicit pole in Eq. (13) 
which leads to a delta function 8{x\—x). Because this pole results into zero gluon momentum 
insertion to the diagrams, we can combine these contributions with the second term in 
Eq. (14) as we mentioned above. Therefore, we can add them together, 

d 



— \ dxT F (x,x)d l g ±lgj_ 



H(kJ g )j> 



|fc±=o 



H^xP,xP,l g )]i> 



(19) 



From this equation, we find that the contribution from Fig. 1(b) cancels out that from the k± 
expansion on the quark line with momentum u k" , because they have the same color-factor 
but opposite signs. The Fig. 1 (d) and the k±_ expansion on the quark propagator u k — l g are 
also the same but with different color-factor: color- factor for Fig. 1(d) is —1/2N C whereas 
that for Fig. 1(a) is C F . Their total contributions will add up to a color-factor CU/2. The 
same color-factor (7^/2 appears for Fig. 1(c). Thus, the final result for this contribution will 
be proportional to C^/2. By applying the identity of Eq. (17) again, we can re-write this 
part of contribution as 

d 



- — — / dxT F (x,x)d l g± 



Os_Ca 

"2tt 2 

a s Ca 
'2tt 2 

'2tt 2 



H (xP,l g± ) 



dxT F (x, x)d 2 l g ±Ho(xP, l g ±) 



— - — — / dx 



dl 



9-L 



l 2 



1 + Z 2 

l-z 



T F (x,x) , 



(20) 



where H represent the hard scattering part without color factor and we have made use of 
the fact that the hard part H (xP,l g ±) oc l/^j_- Again, the same splitting kernel appears. 



Finally, there is also contribution from Tp(xi,X2), which only comes from the hard pole 
diagram Fig. 1(c). Summing up all contributions, we obtain the scale evolution equation for 
the diagonal part of the quark-gluon correlation function Tp(xi,X2), 





<91n [i 2 



T F (x B} fi 



2tt 



r dx 






j X 





1 + z 2 3_. 



C A (1 + z 

~2~ I F ^ ' x > 



1 + z 2 
1 - z 



Tp(x, x) 



Tp(x, x) + Tp(xz, x) 



(21) 



which is consistent with that in Ref. 



25 



261 ] . The complete evolution equation for Tp(x) 



shall also contain contributions from the three-gluon correlation functions, which have been 
calculated in [26]. 



As we mentioned above, the contributions from Figs, l(b-d) to the evolution of g(x) are 
different from that of Tp(x). For g(x) splitting, we do not take pole contributions from 
these diagrams. For example, we will not have cancelation between diagrams Fig. 1(b) and 
collinear expansion of quark line "fc" of Fig. 1(a). More over, without taking pole there will 
be an additional integral variable in the splitting function, similar to that for the evolution 



of Qt structure function 



231 ] . The At contribution from Figs, l(b-d) can be transformed 



into T F and Tp, or to Gp> and Gp, [17J. Because we do not take a pole for the scattering 
amplitudes, the calculations for these diagrams are straightforward. The partial derivative 
contribution from Fig. 1(a) is similar to that for Tp(x) calculation. This part depends on 
g(x). After adding all these contributions together, we obtain the evolution equation for 



dg(x Bl ii 2 
<91n/i 2 



O'.s 



dxdy f 



2vr 

+ G D (x,y) 
+ G D (x,y) 



x 



x 



B 



+ 



— x) 


c F 


X B 


2x 2 B 


y 


xy 


X B 


x B 


y 


X 



1 + z 1 
x B 

X 



C A l + z 2 



z - 



2 1 



- 1 



C A {x\ + xy)(2x B -x-y) 



(x B -y)(x- y)y 



1 



C A x\ 



xy 



2 (y-x B )y_ 



(22) 



where again z = x B /x, and the definitions of Go and Gd follow that in |l7j. The end-point 
singularity from g(x) with color factor Ca/2 at right hand side of the equation is canceled 
out by that from at the second line. We further notice that we can replace Gn and 
Gp, with Tp and Tp at the right hand side by using the relations between them 15|, [19]. 
However, we still have the g(x) term at the right hand side of equation. Although we can 
re- write g(x) in terms of and Tp [yj], that will not eliminate its dependence completely 
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and the right hand side will depend on Go, Tp and Gp, instead. Therefore, the evolution of 
g{x) depends on three functions: g(x), Gp,(x, y) and Gp>(x, y). This feature is different from 
that for Tp(x), where it only depends on Tp and Tp. It may indicate the nontrivial QCD 
dynamics associated with the evolution of the correlation function g(x). This has also been 
shown in its contribution to the Drell-Yan dilepton azimuthal asymmetry in pp scattering. 
We leave this study in a future publication. 

Since the derivation follows the similar procedure, we skip the technique details and only 
list the final result for the evolution equation of correlation functions Tp (x,x),h(x). For 



Tp \ we have 



8^ Tf (W) - 27 



dx 

x 



C F 

c A 



2z 



1-* + 



+ 25(1 



Tp (x, x) 



+ "TT { -r^TP(xz,x) 



1 - z 



2z 
1-z 



Tp (x, x) 



(23) 



which is consistent with the large transverse momentum Boer-Mulders function h^(x,k±) 
calculated in [l{|. Accordingly, we obtain the evolution equation for h, 

dh(x B ,fJ?) 



<91n/i 2 



a s f dxdy J 
+H D (x,y) 



h(x)8(y — x) 



C h 



2z 



+ 25(1 



C, 



2(x - y - x B ) 



+ 



Ca 2xb(xbx + x B y 



C A 2z 



X 



2 1 

y 2 Y 



(24) 



y 2 (x B -y)(x-y)y 

where the twist-three function Hp>(x, y) has been introduced in the Ref. \v\ - Similar to that 
of g(x), the evolution of h depends on h and Hp,. 

In conclusion, we have derived the scale evolution for the transverse momentum dependent 
quark-gluon correlation functions associated with the four k±-odd TMD quark distributions. 
We have performed our calculations in light-cone gauge with a particular boundary condition 
for the gauge potential, and we have checked that our results do not depend on these 
boundary conditions. Our result on the evolution of Tp{x) confirms recent calculations 25, 



261 ] . The scale evolution for g and h reveals nontrivial QCD dynamics. We hope this will 



stimulate further theoretical studies. 

Meanwhile, we notice that the scale evolution for the general twist-three operators have 
been calculated in the literature 23] . It will be interested to compare the evolution equations 
for the correlation functions studied in this paper with these well-known results. Especially, 
the evolution of the twist-three distribution gx{x) and its contribution to semi-inclusive 



11 



processes deserve further investigations. We will address these issues in the forthcoming 
papers. 
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